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The derivation of the response function of an interferometric gravitational wave detector is a 
paradigmatic calculation in the field of gravitational wave detection. Surprisingly, the standard 
derivation of the response wave detectors makes several unjustifiable assumptions, both conceptual 
and quantitative, regarding the coordinate trajectory and coordinate velocity of the null geodesic the 
light travels along. These errors, which appear to have remained unrecognized for at least 35 years, 
render the "standard" derivation inadequate and misleading as an archetype calculation. Here we 
identify the flaws in the existing derivation and provide, in full detail, a correct derivation of the 
response of a single-bounce Michelson interferometer to gravitational waves, following a procedure 
that will always yield correct results; compare it to the "standard", but incorrect, derivation; show 
on : where the earlier mistakes were made; and identify the general conditions under which the "stan- 

dard" derivation will yield correct results. By a fortuitous set of circumstances, not generally so, the 
final result is the same in the case of Minkowski background spacetime, synchronous coordinates, 
transverse-traceless gauge metric perturbations, and arm mirrors at coordinate rest. 
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Ctf ' PACS numbers: 04.80.Nn, 95.55.Ym, 04.30.Nk 
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The derivation of the response function of an interferometric gravitational wave detector is a paradigmatic calcu- 
^Jy lation in the field of gravitational wave detection. The same general derivation appears, almost without variation, in 
i— i \ reference books [l[ , pedagogical and review articles 0, 0, 13| , general relativity texts @ , proposals for experimental 
facilities and technical reports @, 0], and other highly-cited articles from the primary literature (e.g., (1. M. [Tol Hi}). 
V/ | Derivations of the response come in two flavors. By far the most popular is based on the vanishing proper length 
along a null geodesic separating the ends of each interferometer arm. The second flavor of derivation, whose validity is 
' restricted to wavelengths much larger than the size of the detector, uses the equation of geodesic deviation to calculate 
. the changing separation between the mirrors that bound the detectors arms. 

Surprisingly, both these derivations of the detector response function are conceptually and methodologically in- 
correct. Their flaws appear to have gone undetected for at least 35 years [19(, going back to the original proposals 

■ for constructing gravitational wave detectors of this type [(| [T3|. Perhaps equally surprising, a correct derivation 
of the response, which exposes the several errors made in the classic derivations, also finds that the combination of 

■ errors cancel in the special case of gravitational wave perturbations about flat spacetime in synchronous coordinates 
and transverse-traceless (TT) gauge. In more general circumstances, however, the errors are significant and using 
the standard derivations as a template for similar calculations in these more general circumstances will lead to order 
unity fractional errors. 

In Sec. [IT] we recap the standard derivations of the response of an interferometric gravitational wave detector and 
highlight the conceptual errors that invalidate them. In Sec. Mil we provide, in full detail, a correct derivation of the 
response of a single-bounce Michelson interferometer to gravitational waves. We take special care to justify every step 
in the derivation and highlight where the intermediate results in the standard calculation are quantitatively incorrect. 
Finally, in Sec. II VI we discuss the special circumstances that lead the intermediate errors to cancel when the entire 
calculation is assembled and applied to find the detector response. 
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FIG. 1: Spacetime diagram of the optical paths in a simple Michelson interferometer. The dashed lines lie in a constant time 
hypersurface. The observer measures the difference in the phase of light that has traveled along the two different optical paths 
(A) and (B). In the case of a common, monochromatic light source at the beam splitter this difference is the difference in phase 
of the light source at times t — ta and t — tb- 



II. "CLASSIC" DERIVATIONS OF THE DETECTOR RESPONSE 
A. Physical observable 

Consider a simple Michelson laser interferometer, shown schematically in Fig. [TJ By various devices (cf. one 
can measure the quantity (A(f>(t)) mod (2tt), where A0 is the phase difference 0(a)(4) — 0(b)(4) m the light arriving 
at the interferometer output port along optical paths (A) and (B). The light arriving and recombined at the beam 
splitter along these paths originated at the beam splitter at earlier times t — ta and t — tb , so that 

A0(4) = (A) (t) - ( B)(i) = 0o(4 - ta) - 0o(4 - r B ) (2.1) 

where 0o is the phase of the light leaving the beam splitter as a function of time. If the light source is, e.g., a 
laser with frequency v and we choose coordinates such that 4 is proper time along the beam splitter's world line 
then A(f>(t) = 2~kv(t& — ta). In this way a laser interferometer is sensitive to gravitational wave perturbations in the 
difference in time, measured at the beam splitter, required for light to propagate along the optical paths (B) and (A). 
To determine the interferometer response, then, is to determine the gravitational wave perturbation to the difference 
in these light travel times. 

B. Derivation from vanishing proper length of null paths 

The most popular derivation of the interferometric gravitational wave detector response is based on the vanishing 
proper length along a null path separating the two ends of each interferometer arm. Examples of this derivation are 
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found in reference books JT[ , review articles , pedagogical articles and summer school procee ding s E l H| , proposals 
and technical reports [1, U and other highly cited articles from the primary literature [|[ [tj [l(| Illl.ll2j]. The starting 
point for this derivation is the invariant interval in the perturbed spacetime evaluated in transverse-traceless (TT) 
gauge [l3T ] for a null path connecting the beam splitter and an end mirror: 

= g^dx^dx" = -dt 2 + 5 lJ dx i dx j + h i3 dx l dx J . (2.2) 

This expression is rewritten as 



dt = I dX \l dX dX (ftj + ft<jO 



U P1 d\ dX / 

where the perturbation hij is taken to be small and the integration path is along the null geodesic. The integral on 
the left is At, the light travel time along the null geodesic path between the beam splitter and an end mirror. To 
evaluate the right-hand side integral one introduces the unperturbed spatial trajectory of an electromagnetic phase 
front as it propagates alonga null-ray from the beam splitter to the end mirror at z J (cf. e.g^ [l|, Eq. 2. 2. 151, El, 
Eq. 1-3], [| E q s - 2.3-2.4], J| Eqs. 2.16-2.18], @, pg. 2], @, pgs. III-1-III-3] [H Eqs. 35-39], [| Eqs. A2-A4], [lfl, 
Eqs. 4-5]): 

t = A, x j = \( j (2.4a) 
1 = ?t k 6 jk . (2.4b) 

With this substitution the integral on the right becomes 

£d\ (l + ^wfo^XV) = L+ l - j\\h kl {t,x 3 )Q k Q l . (2.5) 

When the gravitational wave period J" 1 is much less than L (i.e., JL <C 1) the perturbation hij is constant along the 
path and the second integral simplifies to LhijQC? /2, where Q is the unit vector coordinate direction from the beam 
splitter to the end mirror. The result is 



At = L+^hijCC j . (2.6) 



What is observed at the beam splitter is proportional to the difference in the light travel time, from the beam splitter 
to an end mirror and back, along each arm, or 

6(At) /.//,, (CO' (2.7) 

where is the unit vector along the second interferometer arm, also of length L, and one assumes a single round- 
trip of the light along each arm. The generalization of this result to a multibounce delay line, Fabrey-Perot cavity 
interferometer, and/or a light storage time greater than the wave period is straightforward. 

A careful reader may question the evaluation of the right-hand side integral in Eq. (|2.3p along the unperturbed 
null-geodesic trajectory described in Eq. (|2.4|) leading to the right-hand side of Eq. (|2.5|) . In particular, use of the 
unperturbed spatial trajectory of the null-ray ignores 0{h) corrections to the coordinate trajectory dx 3 /dt and the 
coordinate x 3 . Why these two O(h) corrections should be ignored in the evaluation of the 0(1) contribution to the 
integrand in Eq. (|2.3p is not apparent and never justified. That the spacetime trajectory is geodesic is not relevant 
here: the integral in question is not over the spacetime path and, in any event, the coordinate trajectories of the null 
geodesic connecting the beam splitter and end mirror world lines in the perturbed and unperturbed spacetimes have 
different end-points — indeed, it is exactly the difference L — r that we desire to calculate. 

In fact, the standard derivation is flawed on exactly this point. Naively following this "classic" derivation in 
circumstances where it may, on the basis of all standard presentations that have been examined, be thought to apply 
(e.g., availability of TT gauge and coordinate stationary mirrors, which are the usually remarked- upon circumstances) 
will, in general, yield incorrect results. As discussed in Sec. IIV1 for perturbations of a restricted set of spacetimes in 
special coordinates the errors resulting from the use of the unperturbed trajectory cancel; however, the reasons for 
this cancellation are nontrivial and have not previously been remarked-up on or examined. 
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C. Derivation from geodesic deviation 

A second type of derivation, based on the Eq. of geodesic deviation [14 1. can be found in the primary literature 
description of interferometric detectors [HJ [1, Eqs 2.4-2.8], textbooks [3, Eqs. 9.24-9.28], and review articles [J, 
Eqs. 2.12-2.14]. In this derivation one identifies a smooth, one-parameter family of geodesies 7a(^): i-e., for each 
A 6 I. 7a (t) is a geodesic with affine parameter t. The vector field (^(X) = d/dX connects neighboring geodesies 
(e.g., the beam splitter world line and the end mirror world line) and also satisfies the Eq. of geodesic deviation, 

= -R-^CU 8 = -irrC 7 |^ 7 + 0(h 2 ). (2.8) 

The geodesies 7a are identified to be the coordinate fixed world lines of the interferometer end mirror and beam 
splitter, and the variations in the length of £ are equated to variations in the light travel time between the beam 
splitter and the end mirror, leading to 

|C , U 1/2 = i(l + ^-<*C J '), (2-9) 

valid in the limit that fL <C 1 , where / is the gravitational wave frequency. 
In this argument one makes the following assumptions: 

• The spatial projection of the null-ray separating the nearby geodesies is also the geodesic deviation vector 

• The coordinate speed of light propagation is unity so that the spatial coordinate distance traveled along £ M (i.e., 
L) can be equated to the coordinate propagation time. 

Both of these assumptions are false at order O(h). Incorrectly assuming that is the spatial projection of the 
null-ray tangent leads to O(h) errors in the spatial distance traveled by the null-ray. Incorrectly assuming that the 
coordinate velocity of light is unity also makes errors of order 0(h). 

In a pe dag ogical article on the importance of relating physical observables only to gauge invariant quantities, 
Garfinkle [15[ finds, in the limit of slowly varying fields, a second-order differential Eq. for the phase difference A<fi 
in terms of the Riemann tensor. This Eq. has the form of the geodesic deviation Eq. in this limit; however, derived 
in the way shown it avoids the questions that arise when the geodesic deviation Eq. is invoked from the outset and 
without argument or justification, and can be correctly generalized to more general circumstances. 

III. THE RESPONSE OF A FREE-MASS, ONE-PASS INTERFEROMETRIC GRAVITATIONAL WAVE 

DETECTOR 

A. Introduction 

In this section we derive the response of a one-pass interferometric gravitational wave detector. While the deriva- 
tion's details will be specialized for the case of free masses and Minkowski background spacetime we take special 
care to indicate where the different specializations are made and how the calculation would proceed in the more gen- 
eral case. The presentation that follows thus serves as a template for more general applications that require similar 
calculations. 

As described in Sec. Ill Al the observable is the difference in phase A<fi of the electromagnetic wavefronts inbound 
at the beam splitter. As shown schematically in Figure [TJ the phase of each inbound wavefront is equal to the phase 
of the electromagnetic wave outbound from the beam splitter at an earlier time. The principal goal in calculating 
the response of an interferometer is thus calculating the elapsed proper time at the beam splitter between when an 
electromagnetic wavefront is inbound at the beam splitter and when it was earlier outbound from the beam splitter. 
If the beam splitter's spacetime trajectory is z /i (r), where r is the proper time measured along the beam splitter 
worldline, the elapsed coordinate time At is readily converted to the elapsed proper time At either by solving the 
implicit equations 

t.\t-At) = I -At, (3Ja) 

or by evaluating the integral 



t-At 



IT 



At = I dt [ — ) . (3.1b) 
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Henceforth we concentrate on calculating the coordinate time difference At separating the emission of a wavefront 
and its subsequent return to the beam-spliter. We do so in two steps. In the first step we calculate the coordinate 
time required for an electromagnetic wavefront to propagate from the end mirror to the beam splitter, arriving at 
time t; in the second step, we calculate the coordinate time required for the wavefront to propagate from the beam 
splitter to the end mirror, arriving at time t' . From these the coordinate time of the round-trip can be calculated for 
each arm and, as described above, the observable phase difference Ac/)(t). 



B. Inbound phase-front propagation 



Focus attention on an interferometer arm consisting of a beam splitter and an end mirror. An electromagnetic 
wavefront arrives at the beam splitter from the end mirror at time coordinate time t. Define Ai(t) as the elapsed 
coordinate time between the event of the wavefront leaving the end mirror and arriving at the beam splitter at time 
t; i.e., if the wavefront arrives at the beam splitter at time t it left the end mirror at coordinate time t — Aj(t). 



Trajectory 



The trajectory x M (A) of an electromagnetic wavefront can be evaluated from the geodesic equation: 

d 2 x» „ dx a dvP 
dX 2 a P dX dX 

where A is an affine parameter along the geodesic and the condition that the geodesic be null is 

dx» dx v 
° = ^dA IX- 



(3.2a) 



(3.2b) 



We introduce the gravitational wave perturbation by writing the metric as the sum of a background gffi and a 
gravitational wave perturbation eh^ u 



9^ = 9 { $+eh liV + 0{e 2 ), 



(3.3) 



where we are explicit about the perturbation expansion parameter e. 

The first step in our calculation is to solve equations (|3.2[) . subject to the constraint in Eq. (|3.2b[) , perturbatively 
about e = 0; i.e., writing the trajectory x M (A; e) as 



(3.4a) 



we solve equations (|3.2[) order-by-order in e. 

For the particular case of interest the background spacetime is Minkowski (i.e., gffl = rj^), the beam splitter and 
the end mirror are free and, in the limit e — > 0, at relative rest. We take advantage of this circumstance to express 
the gravitational wave perturbation in the corresponding TT gauge. In this spacetime and with these coordinates the 
connection coefficients arc given by 



and the geodesic equations become 



d 2 t e dx 1 dx m dhi 
= — + 



dX 2 2 dX dX dt 



0(e 2 



= 



d 2 x> 
dX 2 



dxP 
~dX 



sjl dt dhi p dx™ dh[ p j. jm 1 dx 1 dh lp 



where we have written t for x*. 
At order e° we have 



= 
= 



ePto 
dX 2 

dX 2 



dX dt 



dX 



dX dx r < 



2 dX dx" 



0(e 2 



^0 



dXg j 

dA = ffon ° 



to = cr (A - Ao) 

xj = (A - Xo)a Q n 3 + y J 



(3.5) 

(3.6a) 
(3.6b) 

(3.7a) 
(3.7b) 
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where Ao, &0: n 3 and y 3 are all constants of the integration, with Ao and ao corresponding to the arbitrary zero and 
scale of the affine parameter A, and the n J and y? corresponding to the wavefront propagation direction and location 
at A = Ao- The condition that the trajectory is a null-geodesic is enforced by Eq. (|3.2b[) , which becomes at O(e ) 

= -l + n l n^S lm 

Now turn to the order e 1 trajectory. The order e 1 corrections to the geodesic trajectory satisfy 

d 2 tx 



(3.7c) 



dX 2 
d 2 x{ 

"dX 2 



dh h 



^6 i m Un lm 

" ' bit ' 1 1 



dX 



j dto dhi p , ^ tixo™ dhi p 1 XJ . m dx 0/i, p 



dX dt 



dX dx m 2 dX dx 



(3.8) 
(3.9) 



where | indicates that the partial derivatives of h\ m are evaluated along the Oth-order trajectory (tcxj). For plane 
wave him we have 



him(t, x j ) = hi m (u(t, x)) 
u(t, x) — t — kjX 3 



where k l is a unit vector in the unperturbed geometry 
Consequently, 



= -l + k l k m 5 lm 

dhim _ dh lm du _ dh ir 

dt 
dhi„ 



du dt du 
d\\ lm du 



dht 



du 



dxi du dxi 
Writing for the metric perturbation along the wavefront 

H lm (X) = h lm (u(t(A) |X (A))) = h lm (t(A),x*(A)) 

we have 



dhi 



dHu 



du dX 

Invoking the O(e ) trajectory (eqs. 13 . T[) we find: 



dX 



u(t(A),x>'(A)) 



(J ,,,, ,„,, du dt du dx? / , ., 
Tx u(t(X)^(X)) = -- + — — = a (l - kj ni 



O(e). 



Thus we may write 



dh 



in I 



1 



dHi„ 



dt 



cr [1 - kjTll 

dhi m kj 
dxi ~ 



dX 

dHi, n 



do (1 — fcp^o) dX 



0(e) 
F0(e). 



For plane wave hi m we thus have 
dt\ 

ti(A) = a 
dx{ 

lix =a ° 



1 rr 

7 h 



1 n l n m 

a x (A - Ao) - - , H ; . m (Ao,A,to,x J ) 



n 3 



2 1 ^p^Q 



1 nl 5i m k 



(A - A )n? - + - 1 °_ ^ j ngW ip (Ao,A,to,x J ) 



(3.10a) 
(3.10b) 

(3.10c) 

(3.11a) 
(3.11b) 

(3.11c) 
(3. lid) 

(3.11e) 

(3.11f) 
(3-llg) 

(3.12a) 
(3.12b) 

(3.12c) 
(3.12d) 
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where y{, n\ and o\ are constants of the integration and we have defined the functional TL of the trajectory {t(A), x? (A)} 



as 



H im (A ,A,t,x J ) = / dX'H lm (X') = dX'h lm (t(X'),x. j (X')). (3.12e) 



Finally, the condition that the trajectory be a null-geodesic becomes 

dt dti dx l rfx™ dxj, dx™ 

^ = 77T77T + TTT^TT /m + TTT^TT im ( to ' o) 



g?A c?A c?A dX dX dX 

(n l n?5 lm - at) al (3.13a) 



2. Renormalization 

It will in general be the case that eti and exj involve terms that are unbounded in A. For example, in the calculation 
presented here the terms in question are those linear in (A — Ao). Owing to this unbounded growth the perturbation 
expansion breaks down when (A — Ao)e is no longer small. These terms are generally evidence that "constants of the 
integration" introduced in solving for to and Xg should no longer be considered constants at 0(e). The next step in 
determining the null-geodesic trajectory of an electromagnetic wavefront is to regularize any divergent terms in the 
perturbation expansion for t and x J . This is generally done using renormalization group methods [l6l ll7T ] . In our 
particular case (Minkowski space and TT gauge perturbation) the renormalization is particularly simple: introducing 
the renormalized a, n J , and y J , 



a n J 



cr = cr(l - eo-i) 
X - enj) 



Vo 



(3.14a) 
(3.14b) 
(3.14c) 



the renormalized perturbed null-geodesic trajectory becomes 



dt 

dX 



1 - 



n'n r 



2 1 - fcjiV 



H lm (X) 



t = a 



(X - A ) - 



n n 



2 1 - fc,n9 



+ 0(e 2 ) 
7~tlm(Xo, A, to,x ) 



dx? 
~dX 



n 3 - e n'(5 J 



1 n'n r 



2 1 



S 3p k p H lm (X) 



x J = y 



(A - A ) -v> - e n fc ^' + - 



1 



k 1 

n n 



2 1 



k q ni 



V m k m W H (A 0) A,t 0) xj) 



+ C(e 2 ) 



and the condition that the geodesic be null becomes 

O = -l + n'n m «S ;m + 0(e 2 ). 
With this renormalization the perturbation of the null-geodesic trajectory is seen to always remain small. 



(3.15a) 
(3.15b) 
(3.15c) 
(3.15d) 

(3.15e) 



3. Coordinate velocity along the null trajectory 

Equations (|3.15p give the geodesic null-ray trajectory followed by the electromagnetic wavefront as it propagates 
between the end mirror and the beam splitter. As suggested in Sec. [II] it is not the unperturbed trajectory (A, A^- 7 ); 
rather, it involves 0(e 1 ) corrections in both magnitude and direction when compared to that coordinate trajectory. 
Additionally, the coordinate velocity, 

= l-en l n m H lm + 0(e 2 ), (3.16) 
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also involves 0(e) corrections. There is no a priori reason to expect that the contribution 



. dx 1 dx 3 

^dF"* (3 ' 17) 

to the integral presented in Eq. (|2.3[) should evaluate to L + 0(e 2 ), as is glibly assumed in the classic derivations when 
the 0(e) corrections to the spatial projection of the null-ray, and the 0(e) corrections to the coordinate velocity, are 
ignored without comment. 

4- end mirror and beam splitter worldlines 

The null-ray geodesic trajectory described above intersects the worldlines of the end mirror and beam splitter at 
coordinate times U and tf. Our immediate goal is to determine Ai(t/) = tf — fy. To achieve this goal we must know 
the beam splitter and end mirror worldlines to 0(e 2 ). We can solve for these worldlines using the same methodology 
as we used in Sec. IIII BTHIII B 2l to find the null-geodesic trajectory of the electromagnetic wavefronts. Now, however, 
the equations to be solved for the end mirror or beam splitter worldline z are 

„ d 2 z^ u dz a dz? 

where a M is the four-acceleration acting on the end mirror or beam splitter and boundary conditions are specified on 
z M and dz^ /dr at some initial moment of time. 

In the present case — Minkowski background, TT gauge perturbation, beam splitter and end mirror free and 
initially at coordinate rest — the four-acceleration vanishes and dz* 1 /dX = 5^ at the initial moment of time. The 
resulting (geodesic) equations for the beam splitter and end mirror worldlines have the trivial solution 

z*( T ) = z 3 +0(e 2 ) (3.19) 

where z 3 is the beam splitter or end mirror location at the initial moment of time. 

5. Coordinate light travel-time 

Given the 0(e 2 ) worldlines of the electromagnetic wavefront (t,x J ), the beam splitter Zg S (t) and end mirror z 3 EM (t) 
in the perturbed space time we can calculate the coordinate time interval A, (t) between the events corresponding to 
the intersection of the null-geodesic (t, x 3 ) with the two worldlines z 3 BS and z^ M by solving the implicit equations 

t(A/)=i (3.20a) 

4s(*) = ( 3 - 20b ) 

4 M (t(A i ))=x^(A i ). (3.20c) 

The interval A,(t) is equal to t — t(A-j). The eight equations (|3.20|) and (13. 15c|) determine the two constants A/, \i 
and the six constant y J , n 3 in Eq. (I3.15d|) for x 3 (A). 

In our particular case Zgg and are constant to (9(e^). Without loss of generality we can orient the coordinates 
so that the beam splitter is at the origin and the end mirror is at z J . We introduce constants L and Q 3 , 

C 3 - | (3.21a) 

= -1 + C'C^Jm, (3.21b) 

which are independent of e. In the absence of the perturbation L is the physical separation, at constant t, between 
the beam splitter and the end mirror. Finally, we choose the constant affine parameter scale factor a [cf. Equation 
13. 14a] to be unity. 

Consider a null-geodesic emerging from the end mirror at Ao = and traveling in the direction 

dr 3 ( 1 n'n™ 1 \ 

' -(0) = n'-e n^ m + -- 5™k p ) H lm (0) + <D{e 2 ). (3.22) 



dX w V 2 1 - ki nJ 
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If the wavefront is to arrive at the beam splitter at affine parameter A then n° must satisfy 

/ 1 n l n m \ 

= 2? + An-' - e (n l P m + - — nq k P P p j H im {0, A,«d,iJ) + 0{e 2 ), (3.23a) 

When e = 0, n- 7 is equal to —C? and A is equal to L. As e is increased n J (e) and A(e) change. Noting, however, that 
n- 7 is always — £ J + 0(e) we have 

= + An-' + e [C l P m - ^ k p S j Pj Hi m {0, L, to, x J ) + 0{e 2 ). (3.23b) 

Contracting this expression with Sjin 1 and SjiQ we find 

= LC l n m 5 lm + A - e (l - I J^-^ j ( l C n H lm (0, L, t , xj) + 0{e 2 ) (3.24a) 

1 

2 1 + k q (i 
which together yield 

X = L + e(l- I yT^ ) C ' rWim( °' L ' io ' X ° o) + ° {e2) - (3 ' 25) 

Combined with our earlier expression for t(A) [cf. Equation 13. 15b| we find that the light reaching the beam splitter 
at coordinate time t left the end mirror at coordinate time Aj(t) earlier, where 

Ai(t) = L + ^( l ( m H lm (t - L,t,t ,xl) + 0(e 2 ) (3.26a) 

t (A) = A (3.26b) 

x^(A) = ( J (t - A). (3.26c) 



= L + \5 lm ( l n m + e ( 1 - \ - kp f ' ) ('^(O, L, to, ij) + 0(e 2 ) (3.24b) 



6. Propagation path and angle- of- arrival scintillation 

In the classic derivation the light ray is assumed to propagate along the coordinate straight line path between the 
beam splitter and the end mirror. From Eqs. (|3.15d|l and (j3.23al) we see that the spatial coordinate followed by the 
light ray emerging from the end mirror at affine parameter A = is the "wavy" line 

x*(A) = z J + AnJ + e (qH^ 1 - \ C -^- 8^k m ) H kl (0, A,t ,x J ) + 0{e 2 ) (3.27a) 
V 2 1 + k q (i J 

where 

n j = _ c , + e (pp _ cv) ^ + 1 C l n lm (o,L,to,4) + Q{e2) _ (3 27b) 

A corollary is that the arrival direction of the inbound null ray is time dependent; i.e., the passing gravitational wave 
causes the distant source of light to undergo angle-of-arrival "twinkling" . 



C. Outbound light propagation time 

In the same way that we obtained the coordinate time A;(t) required for a null-ray to propagate from the end 
mirror and arrive at the beam splitter at time t, we can find the coordinate time A (t') required for a null-ray to 
propagate from the beam splitter and arrive at the end mirror at t' . An electromagnetic wavefront arriving at the 
end mirror at time t left the beam splitter at time t — A (t), where 

A (t) =L+ e -C l C m U lm {t-L,t 1 t Q X) + 0{e 2 ) (3.28a) 

and the path along which we calculate Tii m is now 

t (A) = A (3.28b) 
*o( A ) =C j (X + L-t). (3.28c) 

for t - L < A < t. 
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D. Round-trip light travel time 

Light arriving at the beam splitter at time t will have left the beam splitter at an earlier time t — Art (t) , where 



A RT (t) = At(t) + A Q (i - Ai(t)) 

= A i (t) + A (t-L) + 0(e 2 ). 

Consider the monochromatic gravitational plane wave 



(3.29) 



(3.30) 



where u = t — kjX J , kj is the gravitational wave propagation direction and ei m (kj) is the wave's polarization vector. 
The inbound (to the beam splitter) light propagation time is 



= ei ro (Aj)K / dX' exp [2mf (A' + fc^A' - t))] 



= e Jm sine (nfL(l + kj( j )) 
The outbound light propagation time is 



,2iri/ t-L+L(l-fc 3 C J )/2j 



Ao(i -L)=L- £ -Q l rC<fl{t -2L,t- L, kj,C j ) + 0(t 2 



t-L 



Am (* - = e «m / dX ' ex P [ 27rl / ( A ' - k lC 3 ( X> + 2L ~ *))] 



isinc (7r/i(l-%C J ))^ 



= 27ri/[t-L-L(l + fe 3 C 3 )/2j 



e lm(kj)- 



(3.31a) 



(3.31b) 



(3.32a) 



(3.32b) 



The gravitational wave correction to the round-trip light propagation time along arm Q 3 for light arriving at the beam 
splitter at time t is thus 



A C (*)-2L=- £U(t)+C%(t-L) Creimik,). 



*(°), 



(3.33) 



E. Interferometer response 

Now consider the two arms of an interferometer, £■? and ^ k . The response of the interferometer is defined as half 
the difference in round trip light travel times for light arriving at the beam splitter along the two arms: 



R(t,C,Z k ,e lm ,h) 



A c (t)-A 4 (t) 



(3.34a) 



where Af and A^ are given by Eq. (|3.33|) . When the round trip light travel time in the detector arms is short compared 
to the period of the gravitational waves — i.e., fL C 1 — then we are in the small antenna limit. In this limit the 
response becomes 



RsA{t,( J ,V k ,eim,ki) 



(3.34b) 



These expressions are exactly those found in the "classic" derivation. 

Even though neither of the two most important intermediate steps in the calculation of the response agree with 
the "classic" derivation — i.e., the spatial projection of the geodesic null-ray between the two mirrors is not the 
unperturbed trajectory, and the distance traveled along this path is not the unperturbed distance — the combined 
effect of these two errors cancels at 0(e) leading to the same result as the corrected calculation described above. The 
reason for this is the subject of the next section. 
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IV. DISCUSSION 
A. Introduction 



It might seem a remarkable coincidence that 

• The 0(e) correction to the spatial projection of the geodesic path and 

• The 0(e) correction to the spatial coordinate velocity along the geodesic path 

exactly cancel. To understand the conditions under which this cancellation occurs and the classic derivation yields 
correct results we walk through the classic derivation in its most general form, filling-in the implicit, but unremarked, 
requirements necessary for its correctness. 



B. Spatial path predicts coordinate-time 



The first significant implicit assumption in the classic derivation is that the coordinate-time required to traverse a 
null trajectory (to, x ) in the unperturbed spacetime can be expressed as a function of the spatial path Xg alone. Thus, 
the classic derivation begins by expressing infinitesimal elements of a null path as a relation between coordinate-time 
and coordinate space displacements. Expressing the general line-element as 



ds 2 = - (a 2 - f3jpkj jk ) dt 2 + 2f3 j dx j dt + - , : d.r' d.r' . 



(4.1a) 



where the lapse a, shift /3j and spatial metric 7^ are functions of (t, x 3 ), the relation between coordinate-time and 
space displacements for null paths is 



dt = 



a* 



a 2 



-dx 3 



1 



where 7 j7c and f3 2 are defined by 



51 



7jk 



Y 2 -(3 2 



dx 3 dx k . 



Writing the metric coefficients a, /3j and jjk as the sum of a background and a perturbation, 

a = a (o) + 
= 0(o)j + e/3(i)j 
Ijk = 7(0) jk + £7(i) jk, 

the relation between coordinate-time and space displacements for null paths may be written 



dt = 



/3( ) j dx 3 



1 



(3(o)j(3(Q)k\ dx 3 dx k 
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at 



P 2 o 



dX dX 



B, 



dx 3 
~dX 



Gjk 



dx 3 dx 
~d~X~dX 



f3(o)iP(o)m \ dx 1 dx m 



7(0)ir, 



PI 



dX dX 



0(e 2 



dX 



where the Bj and Gjk are linear in the «(i), /3(i)- and 7(i)jfe, and 



P(o) - ^(o)j^(o)fe7(o) 



5! = 



. Ik 



7(o)jz7( )- 



(4.1b) 



(4.1c) 
(4.1d) 



(4.2a) 
(4.2b) 
(4.2c) 



(4.2d) 



(4.2e) 
(4.2f) 



In the classic derivation one implicitly assumes that Eq. (|4.2d|) together with a spatial trajectory x J (A) allows us to 
calculate the elapsed coordinate-time required for a null (electromagnetic) wavefront to travel from spatial coordinate 
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x J (Ai) tox J (A/) along x J : i.e., 

T(x>,\ i ,\ f )=T (x* , Aj, A/) + eTi(x- J , A i; A/) 

/3( )j dx 3 1 



7i(x>') 



P(o)jP(o)k \ dxJ dx k 
7(o)jfc + -2 



dx- 7 dx J dx fe 



\ 



7(o)in 



P(o)iP{o)m \ dx l dx m 
a (o) - ^(o) / dA dA 



dA dA 



dA 



f/A 



(4.3a) 
(4.3b) 

(4.3c) 



For this to be true the 0(1) terms on the right-hand side of Eq. (|4.2d[) must be independent of time; otherwise it 
would also be necessary to know t(X) (and, thus, t(Xi) and t(A/)). 



C. Spatial path perturbations are irrelevant 

The second significant implicit assumption required for the validity of the classic derivation is that the perturbation 
to the null ray's spatial path must be irrelevant to the computation of the coordinate-time required to traverse the 
path. This condition is subtle: it is not the same as requiring that the spatial path be geodesic: geodesic paths have 
extremal path length. It is also not a statement that can be framed in a gauge invariant way on a spatial hypersurfacc: 
it is about the extremality, with regard to spatial path perturbations, of the coordinate time required to traverse the 
path. 

To see how this requirement arises, write the trajectory traveled by the null ray as a function of the perturbation 
parameter e: 

x^A)=x^ o (A)+ e xj(A) + 0( e 2 ). (4.4) 

The elapsed coordinate-time required for the null ray to travel from x J (Ai) to x 3 (A/) along x 3 is 

T(y? ) = T (xl + exj ) + eT a (x J + ex{ ) (4.5) 
= T Q (xi + ocj) + e2i(xg) + 0{e 2 ) (4.6) 

where we have suppressed the designation of the end-points Ai and A/ . In the classic derivation one asserts 

T(J)=T(4)+0(e) 2 , (4.7) 

or, equivalently, 

To^ o + e xj)=To(x£) + 0( e 2 ). (4.8) 

Eq. (|4.8[) is the condition satisfied by a path Xg when it is an extremum of Tq, which is, in general not a geodesic or 
related to a geodesic of either the spacetime or the constant coordinate-time spatial hypersurface. 

Under what conditions can we expect that 7q have an extremum? From Eq. (|4.3b[) we see that it is necessary that 
the /3(o)j should vanish. Additionally, the 0(1) terms on the right-hand side of Eq. (|4.2d[) must be time independent: 
i.e., the ratio 7(o)jfe/ a (o) m ust be independent of coordinate-time. Finally, the end-points of the perturbation - 

ex{(Ai) and ex{(A/) — should be no more than C(e 2 ): i.e., the world lines of the beam splitter and the end mirror 
should be coordinate stationary. 



D. Conclusion 



The classic derivation of the response of an interferometric detector makes a number of significant, unarticulated 
and generally unjustifiable assumptions regarding the background spacetime, the background spacetime coordinate 
system, and the perturbation coordinate gauge that significantly limit its applicability and render it unsuitable as 
an archetype for understanding or for calculation in more general circumstances. For the classic derivation to yield 
correct results 
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• The background spacetime metric must be expressible in a coordinate system with vanishing shift and where 
the ratio of spatial metric to squared lapse is coordinate-time independent; 

• In this background coordinate system the beam splitter and end mirror world lines must be spatial coordinate 
stationary (but not necessarily geodesic); 

• The gravitational wave perturbation must be expressible in a gauge such that the perturbation to the beam 
splitter and end mirror world line coordinates are 0(e 2 ); and 

• The spatial projection of the background spacetime null-geodesic followed by the electromagnetic wavefront 
between the beam splitter and the end mirror must be an extrcmum of the coordinate-time required to traverse 
the null ray [cf. Equation 14. 3b| . 

These conditions hold for flat background spacetimes in Minkowski coordinates, with the beam splitter and end mirror 
at coordinate rest, and gravitational wave perturbation expressed in TT gauge, which is the usual case considered. 
The two most commented upon conditions identified in the classic derivation — that the mirrors are free and at 
relative rest, and that the perturbation is in TT gauge — are in fact neither necessary nor sufficient conditions for 
the validity of the classic derivation (though they are consequences of the necessary conditions when applied to the 
physical conditions of this particular problem) . 

The subtlety of these conditions becomes apparent when other circumstances, where the classic derivation might 
be expected to hold, are considered. For example, consider isotropic cosmological spacetimes. In these spacetimes the 
gravitational wave perturbation may be expressed in a TT gauge with respect to an end mirror and beam splitter at 
rest with respect to the cosmological fluid. When the background is expressed in "arc-time" coordinates, 



ds 2 = a(tf 



-drf 



dr 2 



1 — kr 2 



+ r 2 (d9 2 + 



sin 



2 ) 



(4.9) 



all of the conditions described above hold and the standard derivation will provide the correct result; however, if the 
same physical problem is expressed in the usual Robertson- Walker (synchronous) coordinates, 



ds 2 



-dt 2 + a 2 (t) 



dr 2 



1 — kr 



9 2 + sin 2 I 



deb 2 



(4.10) 



then 7jfc/a 2 is not independent of coordinate-time and the classic derivation fails, even though the spatial hypersurfaces 
and the spatial coordinate expression of the mirror and electromagnetic wavefront world lines are all unaltered. 

In summary, the adage "measure twice, cut once" seems particularly apt. It is fortunate that, in this case, the 
unarticulated and apparently unrealized assumptions made in the classic derivations of the interferometer response 
cancel. Nevertheless, it is important that errors like these be recognized and corrected in order that future calculations 
upon which experiment or observation rely do not make similar mistakes in less fortuitous circumstances. 
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